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Abstract 

The interaction of (linearized) gravitation with matter is studied in the causal approach 
up to the second order of perturbation theory. We consider the generic case and prove 
that gravitation is universal in the sense that the existence of the interaction with gravi- 
tation does not put new constraints on the Lagrangian for lower spin fields. We use the 
formalism of quantum off-shell fields which makes our computation more straightforward 
and simpler. 
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1 Introduction 



The general framework of perturbation theory consists in the construction of the chronological 
products such that Bogoliubov axioms are verified pQ, jl], [3], [9]; for every set of Wick mono- 
mials Wi(xi), . . . , W n (x n ) acting in some Fock space generated by the free fields of the model H 
one associates the operator Tw lt ...,w n (xij ■ ■ ■ i x n)] all these expressions are in fact distribution- 
valued operators called chronological products. Sometimes it is convenient to use another 
notation: T(Wi(xi), . . . , W n (x n )). The construction of the chronological products can be done 
recursively according to Epstein-Glaser prescription [4], [5] (which reduces the induction proce- 
dure to a distribution splitting of some distributions with causal support) or according to Stora 
prescription [T7] (which reduces the renormalization procedure to the process of extension of 
distributions). These products are not uniquely defined but there are some natural limitation 
on the arbitrariness. If the arbitrariness does not grow with n we have a renormalizable theory. 
An equivalent point of view uses retarded products [21]. 

Gauge theories describe particles of higher spin. Usually such theories are not renormal- 
izable. However, one can save renormalizability using ghost fields. Such theories are defined 
in a Fock space "H with indefinite metric, generated by physical and un-physical fields (called 
ghost fields). One selects the physical states assuming the existence of an operator Q called 
gauge charge which verifies Q 2 = and such that the physical Hilbert space is by definition 
"Hphys = Ker(Q)/Im(Q). The space T-L is endowed with a grading (usually called ghost number) 
and by construction the gauge charge is raising the ghost number of a state. Moreover, the 
space of Wick monomials in 7-L is also endowed with a grading which follows by assigning a 
ghost number to every one of the free fields generating "H. The graded commutator d,Q of the 
gauge charge with any operator A of fixed ghost number 

d Q A=[Q,A] (1.1) 

is raising the ghost number by a unit. It means that g?q is a co-chain operator in the space of 
Wick polynomials. From now on [■, ■] denotes the graded commutator. 

A gauge theory assumes also that there exists a Wick polynomial of null ghost number T(x) 
called the interaction Lagrangian such that 

[Q,T]=td,T» (1.2) 

for some other Wick polynomials T M . This relation means that the expression T leaves invariant 
the physical states, at least in the adiabatic limit. In all known models one finds out that there 
exists a chain of Wick polynomials T M , T 7 ^, T Myp , . . . such that: 

[Q, T] = id^T^, [Q, 7"*] = id u T^, [Q, T""] = id p T^», . . . (1.3) 

In all cases T^ v , T MI/p , . . . are completely antisymmetric in all indices; it follows that the chain 
of relation stops at the step 4 (if we work in four dimensions). We can also use a compact 
notation T 1 where J is a collection of indices / = [vi,...,u p ] (p = 0, 1, . . . , ) and the brackets 
emphasize the complete antisymmetry in these indices. All these polynomials have the same 
canonical dimension 

ujiT 1 ) = u , \/I (1.4) 
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and because the ghost number of T = T is supposed null, then we also have: 

ghiT 1 ) = (1.5) 
One can write compactly the relations (II. 3p as follows: 

d Q T J =i9/". (1.6) 

For concrete models the equations (11.31) can stop earlier: for instance in the case of gravity 
rpiivpo — o. In [11] - [TJ] we have determined the most general solution of (11 .31) for the case of 
massless gravity interacting with Yang-Mills, scalar and Dirac fields; in [15] we have considered 
the case of massive gravity in interaction with matter. 

We can construct the chronological products 

T h >~' In {x x , ...,x n ) = T(T h ( Xl ), T In {x n )) 

according to the recursive procedure. We say that the theory is gauge invariant in all orders of 
the perturbation theory if the following set of identities generalizing (11.61) : 

n f) 

d Q T h >~> In = i y^(-l) Sl T h '-' Illl '-' In (1.7) 

are true for all n G N and all Ji, . . . , I n . Here we have defined 

i-i 

= ( L8 ) 

(see also [3]). In particular, the case I\ — . . . — I n — it is sufficient for the gauge invariance 
of the scattering matrix, at least in the adiabatic limit. 

Such identities can be usually broken by anomalies i.e. expressions of the type A 11 '"'' 1 ™ which 
are quasi-local and might appear in the right-hand side of the relation (II. 7p . One compute these 
anomalies in lower orders of perturbation theory and imposing their cancellation one obtains 
various restrictions on the expression of the interaction Lagrangian. In this paper we consider 
the interaction between gravitation, Yang-Mills fields and matter (described by scalar and Fermi 
fields). The complete expression of the Lagrangian has been given in a previous papers [UJ - 
[TJ] ■ Here we compute the anomalies of this model in the second order of perturbation theory. 

We find usefull to use the formalism of off-shell fields for the computation of the anomalies. 
This formalism give a systematic way of computing the anomalies. Such a formalism was used 
previously in the literature in the context of classical field theory. We use here a pure quantum 
version. 

In the next Section we remind our definition of free fields. We avoid explicit formulas using 
the reconstruction theorem of Wightmann and Borchers algebras. In Section [3] we recall the 
main results concerning the interaction Lagrangians for various models with higher spin fields 
(Yang-Mills and massless gravitation). In Section H] we introduce the off-shell formalism. Then 
in Section [5] we describe the consequences of the cancellation of the anomalies in the second 
order of the perturbation theory in the most general case. 
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2 Free Fields 



2.1 Free Scalar and Dirac Fields 

We will adopt the description of free quantum fields given by the reconstruction theorem from 
axiomatic field theory [16] , [22] based on Borchers algebras. In this approach one can construct 
a quantum field giving the Wightmann n-points distributions and the statistics. For a free field 
it is sufficient to give the Wightmann 2-points distribution and generate the rest according to 
Wick theorem. We use formal distribution notations for simplicity. 

A. The Real Scalar Field. We start with the most elementary case of a real scalar field. 
The field is $(x) and the Hilbert space is generated by vectors of the type 

$(xi)---$(x n ) ft (2.i) 

where ft is the vacuum vector. By definition, the 2-points distribution is 

< ft, $(x 1 )$(x 2 )ft >= -i D%\x x - x 2 ) (2.2) 

where Dm\x) is the positive frequency part of the Pauli- Jordan causal distribution of mass 
m. We assume that the scalar field is a Bose field and the n-points distributions are generated 
according to Wick theorem: for n odd 

<ft,$(xi)---$(s„,)ft >=0 (2.3) 

and for n even: 

< ft, $(xi) ■ ■ • $(x n )ft >= ^ < ft, $(a; CT( i ) )$(x (j(2) )ft > • ■ ■ < ft, $(x <T( „_i ) )$(x (7(n) )ft >;(2.4) 

(7 

here the sum is over all permutations a of the numbers 1,2, ... ,n. We also postulate that the 
field $ is self-adjoint: 

$ f = $. (2.5) 

Then one can construct the Hilbert space % from vectors of the type (12.11) with the 
scalar product < •, • > reconstructed from the n-points distributions given above and the 
self-adjointness assumption. We first define a sesquilinear form in the Hilbert space between 
two states of the form (12.11) by 

< $(x n ) • • • $(xi)ft, $(x n+1 ) • • • $(x m+n )ft > 
=< ft, $(xi) f ■ ■ • <$>(x n ) ] <$>(x n+1 ) ■ ■ ■ $(s m+n )ft >=< ft, $(xi) ■ • ■ $(x m+n )ft > (2.6) 

and one can prove that is positively defined so it induces a scalar product. 

Then the action of the scalar field on states of the form (12 .ip is defined in an obvious way. 
One can prove that the scalar field so defined verifies the Klein-Gordon equation of mass m 

K rn <$> = (n + m 2 )<S> = 0; □ = d 2 = d ■ d = <9 M (2.7) 
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and the canonical commutations relation: 

[$(zi),$(z 2 )] = -iD m (xi-x 2 ). (2.8) 

Because of this commutation relation the writing of a state from the Hilbert space in the 
form (12. ip is not unique. 

Moreover, one can introduce in the Hilbert space H a unitary (irreducible) representation 
of the Poincare group according to 

U AA $(x)U£ = $(A- 1 -(x-a)) 

U Ka tt = tt (2.9) 

(here A e C\ is a proper orthochronous Lorentz transform and a is a space-time translation). 
One can obtain in an elementary way the action of the operator U^ a on vectors of the type 
(12.11) by commuting the operator with the factors *&(xj) till it hits the vacuum and gives the 
identity. In the same way one can define the space and time parity operators. 

Of course, one can obtain very explicit representations for the scalar field, but they will be 
not needed in the following. We only mention that one can define in the same way the Wick 
(or normal) products : <& n (x) : for any integer n (see [21]). 

For an ensemble of real scalar fields $ Q where a takes values in a index set I we only replace 
(O by 

< Q, $ a (x!)$ 6 (x 2 )n >= -i 5 ab - x 2 ) (2.10) 

and we make a corresponding modification of the formula (12. 4p . A complex scalar will be an 
apropriate combination of two real scalar fields. 

B. The Dirac Field. For fields of spin 1/2 the preceding formalism remains essentially the 
same with some technical modifications. First one needs a set of so-called Dirac matrices; they 
are 4x4 matrices 7 M verifying the Clifford algebra 

{y,7 I, } = 2iT (2-11) 

where if" is the Poincare invariant tensor (with only diagonal elements which we choose 
1,-1,— 1,-1). Such a set of matrices is essentially unique (up to some similarity transfor- 
mation) according to Clifford algebra theory. In four dimension a simplified analysis due to 
Pauli is available (see for instance [18J). We also define 

75 = 7o7i7273, % = \ {I + e 7s), e = ±, 7 M = ^WY* Y Y). (2.12) 

There exists a four- dimensional representation S(A) of the group SL(2, C) (the universal cov- 
ering group of the proper orthochronous Lorentz group) such that 

S(A)- 1 1 S(A)=A(A)- r , (2.13) 

here SL{2, C) 3 A -»■ A(A) G C\ is the covering map of the proper orthochronous Lorentz 
group and we note that the Lorentz group has a natural action on the Dirac matrices because 
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they carry a Lorentz index. We will denote the matrix elements of these matrices by 7^, a, (3 — 
l,-..,4. 

We define the causal (matrix) distribution 

S M (x) = (i ^ dp + M) D M (x) (2.14) 

and the corresponding positive and negative frequency parts are . Then a Dirac field is a 
set fields ip a (x),a = 1, ■ • • ,4 considered as a column (1 x 4) matrix. It is convenient to define 
the 4x1 (line) matrix ip = ip^ 70. Then the 2-point distributions are given by 

< Q,ip a (x l )i!i3{x 2 )^ >= -i S\y{xx -x 2 ) a /3 

< Q, 4! a {xi)tpi3{x 2 )fL >= -i Sm\x 2 - xi) Pa (2.15) 

and zero in the other two remaining cases. 

The n-points distributions are defined by a formula similar to (12. 4p but we introduce in the 
sum over the permutations the signature (—1)'°"' of the corresponding permutation. This will 
ensure the Fermi character of the Dirac field. One can prove that the field ip(x) fulfills the 
Dirac equation of mass M: in matrix notations 

i ^ d^{x) = M%Ij(x) i d^(x) ^ = -M $(x) (2.16) 

and the canonical anti-commutation relation: 

{ipaiXt), i)p{x 2 )} = -i S M (xi - X 2 )ap- (2-17) 

The unitary representation of the (proper and orthochronous) Lorentz group is defined in 
this case by 

UaM^Wa 1 = SiA- 1 )^ MMA)- 1 ■ x) 

U A tt = tt. (2.18) 

One notices that there is no self-adjointness assumption on the Dirac field. There is however 
an spinor analogous of the real scalar field, namely the Majorana field. One can prove the 
existence of the so-called charge conjugation matrix C defined according to 

7f = -C 7m C" 1 (2.19) 

i.e. the similarity transformation connecting the Dirac matrices with their transposed. Then 
one defines the charge conjugate spinor of ip by 

tfj c = C i) T (2.20) 

and a Majorana spinor is a spinor equal to its charge conjugate ip c = ip. 

We close by mentioning that the choice of the statistics is crucial for obtaining a positively 
defined scalar product. Should we have made the opposite choice in the preceding two case, 
the result would have been a Fock space with a sesquilinear form, non-degenerated but not 
positively defined. 
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2.2 Higher Spin Fields 

For fields of higher spin one can use the preceding formalism with one major modification: it 
is necessary to introduce ghosts fields, which are fields with the "wrong statistics" . 

C. The Massless Vector Field. In this case we consider the vector space H of Fock 
type generated (in the sense of Borchers theorem) by the following fields: (v M , u, u) where the 
non-zero 2-point distributions are 

< Vt 1 v ,1 (x l )v u (x 2 )VL >=i r]^ D^ + \xx -x 2 ), 

< Q, u(xi)u(x 2 )tt >= -i D { ^\xi - x 2 ). (2.21) 

We also assume the following self-adjointness properties: 

vl = v„, u ] =u, v) = -u. (2.22) 

When we generate the n-point function according to a formula of the type (12.31) and (12. 4|) 
we assume that the field v M is Bose and the fields u, u are Fermi. When defining the unitary 
representation of the Lorentz group we consider that the first field is vector and the last two 
are scalars. Because of the "wrong" statistics the sesquilinear form defined by a formula of the 
type (12.61) will not be positively defined. Nevertheless, because it is non-degenerated, we can 
prove that we have Klein-Gordon equations of null mass: 

□ w m = o Du = Uu = (2.23) 

and the canonical commutations relation: 

[v^x 1 ),v v {x 2 )}=if D ( Xl -x 2 ) 

{u(xi),u(x 2 )} = —i D (xi - x 2 ) (2.24) 

and all other (anti) commutators are null. 

We can obtain a bona fidcE scalar product introducing the so-called gauge charge i.e. an 
operator Q defined by: 

[Q, v») = i d»u, {Q, u} = 0, {Q, u} = -i 

QQ = 0. (2.25) 

Using these relation one can compute the action of Q on any state generated by a polynomial 
in the fields applied on the vacuum by commuting the operator Q till it hits the vacuum 
and gives zero. However, because of the canonical commutation relations the writing of a 
polynomial state is not unique. One can prove that the operator Q leaves invariant the canonical 
(anti) commutation relations given above and this leads to the consistency of the definition. One 
can prove now that the operator Q squares to zero: 

Q 2 = (2.26) 

and that the factor space Ker(Q) / Ran(Q) is isomorphic to the Fock space particles of zero 
mass and helicity 1 (photons and gluons) [10] . 



6 



D. Massive Vector Fields. In this case we consider the vector space 7i of Fock type 
generated by the following fields: {y^,u,u, $) where the non-zero 2-point distributions are 

< VL,v^{xi)v u {x 2 )VL >= i rj^ D ( +\xi - x 2 ), 

< Q, u(xi)u(x 2 )Q >= —i D^\x\ — x 2 ) 

< Q, $(xx)<l»(x 2 )fi >= -i D£\xi - x 2 ) (2.27) 
with m a positive parameter. We also assume the self-adjointness properties 

vl = Vft, u ] =u, u ] = -u, $ f = $. (2.28) 

When we generate the n-point function according to a formula of the type (12.31) and (12.41) we 
assume that the fields v^, $ are Bose and the fields u,u are Fermi. When defining the unitary 
representation of the Lorentz group we consider that the first field is vector and the last three 
are scalars. As in the massless case, the sesquilinear form defined by a formula of the type ( 12. 6p . 
will not be positively defined. We can prove that we have Klein-Gordon equations of mass m: 

(D + m 2 )v^ = (D + m 2 )u = (D + m 2 )u = (D + m 2 )$ = (2.29) 

and the canonical commutations relation: 

[v»(x 1 ),v v M]=irr D m (xi-x 2 ) 
{u(xi),u(x 2 )} = -i D m (xi -x 2 ) 

[$(£!), $(x 2 )] = -i D m ( Xl - x 2 ) (2.30) 

and all other (anti) commutators are null. 

We can obtain a scalar product introducing the gauge charge Q according to: 

[Q,v»]=id> M u, {Q, M } = 0, {Q,«} = -2 (9/ + m$) [Q,<$>]=imu 

Qtt = 0. (2.31) 

One can prove that the operator Q leaves invariant the canonical (anti) commutation relations 
given above and this leads to the consistency of the definition. Also one can prove now that 
the operator Q squares to zero: 

Q 2 = (2.32) 

and that the factor space Ker(Q) / Ran(Q) is isomorphic to the Fock space particles of positive 
mass and helicity 1 (massive vector fields e.g. and Z particles) [TO] . 

E. The (Massless) Graviton. The Fock space is generated by the set of fields (h pu , u p , u a ) 
where h pu is a symmetric tensor field (with Bose statistics) and vector fields (with 
Fermi statistics). The (non-zero) 2-point functions are by definition: 

< Q, h fll/ (x l )h pa (x 2 )Q >= - % ~ (Tjpp T) va + T] up T}^ - T] pu r] pa ) -Do +) ( x i - x ^), 



< Vt,u p (xi)u u (x 2 )Vt >= i T)p, v Dfr + '(xi - x 2 ) (2.33) 



7 



and we assume also that 



h\ u = V> u \ = u pi fi J = (2-34) 
One can prove the validity of the equations of motion 

□ V = □ u p = Uu a = (2.35) 

and the canonical (anti) commutation relations 

i 

[hp„(xi), h pa (x 2 )] = -- (rjftp r} va + r} vp rj^ - 77^ r} pa ) D (x 1 - x 2 ), 

{u p (x 1 ),u u (x 2 )} = i r]^ D (x 1 - x 2 ) (2.36) 

and all other (anti) commutators are null. One defines the gauge charge according to 

[Q,hp, u } = --(dpU„ + d v Up- rip U dpU p ), {Q,u ll }=0, {Q, Up} = i d u hp U) (2.37) 

proves that is well defined and that the factor space Ker(Q) / Ran(Q) is isomorphic to the Fock 
space particles of zero mass and helicity 2 (gravitons) [11]. In the same way one can treat the 
case of massive gravity. 

F. We can combine the cases C, D and A described above as follows. We take / = /1U/2U/3 
a set of indices and we make the following choices: 

(a) for a G l\ we consider the set of fields (v£,u a ,u a ) of null mass; 

(b) for a G I 2 we consider the fields (v%, u a , u a , $ a ) of positive mass m a ; 

(c) for a 6 J 3 we consider the scalar fields $ a of mass ; they are called Higgs fields. We 
define the mass matrix m according to 

m ab = 5 ab m a , a, b G I 2 U I 3 . (2.38) 

The expressions for the 2-points distributions are easily given if we want that the final Fock 
should be a tensor space of copies of massless vector fields, \I 2 \ copies of massive vectors 
fields and copies of scalar fields. Equivalently, we can consider that we have the fields 
(u^, u a , u a , $ a ) for all indices a G I\ U I 2 U ^3 but we have $ a , a G h and v% = 0, u a = 0, u a = 
0, a G I3. In the most general case we add cases B and E: 

(d) for A G I4 (another index set) we consider the Dirac fields ipA- As above we define the 
mass matrix of the Dirac fields by 

M AB = S AB M A , A,Beh. (2.39) 

(e) the (massless) gravitons. 

The matter is described by the Dirac fields together with the Higgs fields. 
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3 Interactions 



The discussion from the Introduction provides the physical justification for determining the 
cohomology of the operator cIq = [Q, •] (where by [•, •] we always mean the graded commutator) 
induced by Q in the space of Wick polynomials. A polynomial p G V verifying the relation 

d Q P = i dtff (3.1) 
for some polynomials p M is called a relative co cycle for d,Q. The expressions of the type 

p = d Q b + id^, (3.2) 

are relative co-cycles and are called relative co-boundaries. We denote by Zq 1 , Bq 1 and Hq the 
corresponding cohomological spaces. In (13. ip the expressions p M are not unique. It is possible 
to choose them Lorentz covariant. We have a general description of the most general form of 
the interaction of the previous fields which combines results from [10J and [IT] . Summation 
over the dummy indices is used everywhere. For simplicity we do not write the double dots of 
the Wick product notations. 

A. Yang-Mills Fields 

Here we combine cases C, D, A and B. 

Theorem 3.1 Let T be a relative co-cycle in the variables (v%, u a , u a , <& a ), a G ii U J 2 U J 3 and 
ipA-, A G Z4 which is tri-linear in the fields, of canonical dimension oj(T) < 4 and ghost number 
gh(T) = 0. Then: (i) T is (relatively) cohomologous to a non-trivial co-cycle of the form: 

tym = fabc Q V afM V bv + U a V% d^U c 

+fabc(®a (fib V cfi + m b $ a U b U c ) 

+ 1 fabc ®a <^b^c+Yl ® V W + ZZ K ® ^ $ " ( 3 - 3 ) 

e=± e=± 

(ii) The relation dQt ym = i c^ty m is verified by: 

f&n = fabc \U a V bu FZ» ~ l ~U a U b d»u)j + f' abc $ a ^ U c + ^ ^ ® 7 ^ « a (3.4) 

fmj T7ie relation dqt^ m = i d v t^ is verified by 

t& = \fabc u a u b Ff. (3.5) 

and we have c?Q^ m = 0. 

(iv) The expressions given above are self-adjoint iff the constants f abc , f' abc and f" bc are real 
and the matrices t e a , i s e a are self-adjoint. 
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Here we have defined the gauge invariants which are not coboundaries 

= d»v v a - 9X, VaeAU I 2 (3.6) 

and 

r a = d^ a -m a v^ VaG/ 2 (3.7) 
and the various constants are constrained as follows: 

(a) first, we have some limitations on the range of the indices because of the definitions of 
the sets of indices Ij, j = 1,2,3 namely: (al) f a b c = if one of the indices is in ^3; (a2) also 
f'abc = if c G ^3 or one of the indices a and b are from Ii, (a3) the matrices t e a are zero if 
a <E I 3 ; (a4) the matrices s e a are zero if a G I\. 

(b) The constants f a b c are completely antisymmetric 

fabc = f[abc[- (3.8) 

(c) The expressions f' abc are antisymmetric in the indices a and b: 

fabc = ~f'bac (3-9) 

and are connected to f a bc by: 

fabc m c = f' cab m a - f' cba m h . (3.10) 

(d) The (completely symmetric) expressions f" hc = f" abc } verify 

f „ J -k f'abc ( m l - m l) for a,bel 3 ,cel 2 
fabc I —fate m 2 b for a,cG/ 2 ,6G/ 3 . 

(e) The matrices t e a and s £ a are connected by: 

m a s e a = i(Mt e a -t- e M). (3.12) 

B. Gravitation 

Theorem 3.2 Let T be a relative co-cycle in the variables {h pvi u pi u a ) which is tri-linear in 
the fields, of canonical dimension uj(T) < 5 and ghost number gh(T) = 0. Then: (i) T is 
(relatively) cohomologous to a non-trivial co-cycle of the form: 

tgr = «(2 h w d»h uX d p h vX + 4 h up dyft V d^h pX - 4 h pX d»h Vp d^h v X 

+2 h px d p h pX d p h - h w d p h d p h -Au p d v u x d p h vX 
+4 d p u u d u u x h pX + 4 d p u u d x u v h pX - 4 d v u v d p u x h pX ) (3.13) 

(ii) The relation dQt gT = i d p t p r is verified by: 

t£ = k(-2u p d u h pX d p h vX + u p d p h vX d p h vX - -u p d p h d p h 

+4 u p d u h px d p h vX -2u p d p h vX d p h vX + u p d p h d p h 
-4 d p u v d v h pX h pX - 4 d p u v d x W v h pX + 4 d x u p d p h up h vX 
+4 d„u v d p h px h pX - 2 d„u v d p h px h pX - 2 d p u x h pX d p h + d v u v h d»h 
-2 u p d u d p u p u v + 2 u p d p (fu a u»-2u p d x u p d p u u 
+2 u p d x u p d p u x + 2 d p u p d x u» u x -2u p d p u x d p u x ) (3.14) 
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(Hi) The relation d-Qt 1 ^. = i d v t^ is verified by: 

t£ = k[2{-u" d x u p d p h uX + u p d\u^ d p h uX + u p d p u x d u h px + d p u p d x u p h uX ) 

-{ji v) + 4 dV d p u u h pX ]. (3.15) 

(iv) The relation dqt 1 ^ = i d p t^ p is verified by: 

tg p = k[2u x d x u p v? v - u p (d p u x d x u u - d u u x d x u p ) + circular perm.] (3.16) 

and we have dQt^ p = 0. 

(v) The preceding expressions are self-adjoint iff the constant k is real. 

The gauge invariants which are not coboundaries are more complicated to write. Beside Up and 
the antisymmetric combination 



we also have the linear part of the Riemann tensor from which we must subtract all traces |TTj . 
C. Interaction between Yang-Mills Fields and Gravity 

Theorem 3.3 Let T be a relative co-cycle depending on the Yang-Mills and gravitational vari- 
ables, which is tri-linear in the fields, of canonical dimension oj{T) < 5 and ghost number 
gh(T) = 0. Then: (i) T is (relatively) cohomologous to a non-trivial co-cycle of the form: 

tint = fab (4V F v hp - h F apa F pu + 4 Up d v u a - h pu <ft <j>l + m a Up u a <%) 

+f cd (hp V p* e er* A - m * + m * h $ a $^ 

+ (hp U - I f]p U h) (d^ e ® Yleip - ^ c £ <8> Yled^). (3.18) 



2 

(ii) The relation dgtint = i dpt p nt is verified by: 



Ct = fab (u p F p ° F hpry + 4 u p Fr F bup -2u^ v a <f, bu +4u v <ft Vb , 

(1 2 2 

- u p d v $ c d"$ d - u v d^ c d»<f> d - m - + mb u » $ a $ 6 

~\u u [{d^ C' <g> Ylell> ~ $ C' ® Yled P *p) + {fi V)] (3.19) 

and we also have 

d Q t^ = 0. (3.20) 

(ii) The constants f ab are non-zero only for a,b G I\ U I2 an d f' cd are non-zero only for 
c,d G ^3; these two matrices are symmetric and commute with the mass matrix: 

fab (m a - m b ) = 0, f' cd ( m c - m d ) = 0. (3.21) 
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The matrices c e verify the identity: 

M e = c~ e M (3.22) 

(Hi) The preceding expressions are self-adjoint iff the constants f a b and f' cd are real and we 
also have 

( c *)t = _ c -e. (3.23) 



In all three theorems the co-cycles t 1 are non-trivial. 

There are different ways to obtain the preceding results. One can proceed by brute force, 
making an ansatz for the expressions T 1 and solving the identities of the type (jl.6p as it is 
done in [19]. There are some tricks to simplify such a computation. The first one makes an 
ansatz for T and eliminates the most general relative cocycle. Then one computes g?qT and 
writes it as a total divergence plus terms without derivatives on the ghost fields. Another trick 
is to use the so-called descent procedure. We present them in the most elementary case of pure 
Yang-Mills fields (only the set I\ is non-void). We have the general form: 

+9abc v a u bd^u c + 9a b %vy b u c + g^v£d ^u b u c . (3.24) 
Eliminating relative coboundaries we can fix: 

f (1) _ _ /■(!) f (2) _ n (3) _ q (2) _ (2) (3 

Jabc ~ Jbaci J abc ~ U ' i/abc ~ U > iJa b c ~ iJ b a C - ^O.ZOJ 



Then we obtain easily: 
where: 



dgT = iu a T a + total div (3.26) 



V 



+(-/£ + fell + fi x l + til) < 9A< 

e, upa &*v v b d a v cp . (3.27) 

Now the gauge invariance condition (11.21) becomes 

u a T a = dpi? (3.28) 

for some expression t M which has, from power counting arguments, the general form 

tP = u a t» + d^u a t a + 8 v u a C (3-29) 

where the polynomial t% u does not contain terms with the factor jf" . Then the relation (13.281) 
is equivalent to: 

d^-ml t a = T a 

e a + dn a + d v t7 = o 

C = C- (3-30) 



12 



One can obtain easily from this system that 

T a = (n + m 2 a )t a . (3.31) 
Writing a generic form for t a it is easy to prove that in fact: 

T a = 0; (3.32) 

from here we easily obtain the total antisymmetry of the expressions f^ c and also we have 

(2) (3) 

9abc = 0- New one can take f^ bc = if we subtract from T a total divergence. As a result we 
obtain the (unique) solution: 

T = fUlKv^v^ - vy b d,u c ) (3.33) 

which is the expression from the first theorem in this particular case i.e. the first line of (13.31) . 

Now we briefly present the descent method in this case. There are two results which must be 
used repeatedly [10] . First, we have a version of the Poincare lemma valid for Wick monomials 
and then we have a description of the cohomology group Hq of dq in terms of invariants: if T 
is a Wick polynomial verifying dq T = then it is of the form T = dqB + T where T depends 
only on the gauge invariants u a , F£ v . 

By hypothesis we have 

d Q T = i d^. (3.34) 

If we apply dq we obtain d^dq T p = 0. Using Poincare lemma so one finds out some Wick 
polynomials such that 

dqT" = % d v T^ v \ (3.35) 
Continuing in the same way we find T^d suc h that 

dqT lH = i dpT^d; (3.36) 

we also have 

ghiT 1 ) = \I\. (3.37) 

It means that T^d 

is a sum of terms of the type r/^ u a u b d p u c i.e. is a coboundary: 
T^l = dqB^d. We introduce in (ET351) and obtain 

d Q (TM - i d p B^) = 0. (3.38) 

Using the description of the cohomology of dq we can easily find that we have: 

TlH = d Q B [ ^ ] + % d p B^ vp] + T<J H (3.39) 

where the last term depends only on the invariants u a , F£ v i.e. 

Tt ] = \ /go ^ "» F r + \ /[2 ]c ^ «« «6 F cpa ; (3.40) 
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We substitute these expressions in f !3.35|) and obtain 

d Q (T^ - i d v BW - t") = (3.41) 

where: 



f = /$ ]c ("a v bu - l -u a u b 8»u)j - f^ ]c u a v hv F cpa . (3.42) 
If we use again the cohomology of dg we can easily find out that in fact: 

T M = d Q B» + % d v B^ u] + #\ (3.43) 
We substitute this in f l3.34|) and we obtain the restrictions 

f (l) __ f W f (2) __ f (2) ( o U) 

J[ab]c~ J[ac\b-> J[ab]c~ J [ac]b 

so the constants /r„2 c , /Si c are in fact completely antisymmetric and 

d Q (T"-id v B' w -t) = (3.45) 

where 

' '[,,/«■] \ ~ ' »/' '''»' "" ' /, (/ /<"<- J ~ ~ ./[,,/«■] 1 /"'/^ 

The description of the cohomology of g?q leads to 

T = d Q B + i + t. (3.47) 

Finally one proves that the last term from the expression t is a total divergence. 

In the general case one can use both methods also. In the descent method we will have 
more gauge invariants in the description of the cohomology group Hq. 



t = fSL (I *V Vbu Fr + u a < d,u c ) - i /£ cl e, vpa < < F?. (3.46) 
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4 The Off-Shell Formalism 

It is known that in the second order of the perturbation theory some anomalies can appear and 
this is due essentially because the Pauli- Jordan distribution D m verifies Klein-Gordon equation: 

K m D m = (□ + m 2 ) An = 0. (4.1) 

but the associated Feynman distribution verifies 

K m £>£ = (□ + m 2 ) D m = 5(x - y). (4.2) 

Let us describe in detail this point. One computes the second order causal commutator and 
finds out that the tree contribution has the following generic form: 

[T h (x),T l2 (y)f ee = J2 I D m (x-y) A^(x,y) + d a D m (x - y) A^^y) 

m 

+d a d? D m {x - y) A2J«W(x, y) ] (4.3) 

where the sum runs over the various masses from the spectum of the model and the expressions 
A 1 ^ 12 , A I T } l ' l2 ' ,a , Am' l2 '^ a ^ are Wick polynomials. Moreover we have from (11.61) the identity 

d Q [T^),T'%)]=*A [T ^ (4.4) 

which stays true if we take only the tree graphs. Now one can find out the corresponding 
chronological products by simply substituting in the preceding expression the causal distribution 
by the associated Feynman propagator: D m — > i.e. 

T^(x, y)*™ = J2 l D F m {x - y) A^(x, y) + d a D F m {x - y) A^<*(x, y) 

rn 

+d a d p D F m [x - y) A^ {aP \x,y) ] (4.5) 

In this way all Bogoliubov axioms are true (in the second order) but we might break gauge 
invariance i.e. the identity ( ll.Tjl for n = 2 

d Q T h > l2 (x,y) = i —T I ^ l2 (x,y) + (-l)l^_T^(x, y) (4.6) 

might not be true. Indeed, let us consider the simplest case I\ = \p], I2 = and suppose that in 
the chronological product T(T^(x),T(y)) tree we have a term of the type d^D^x — y) A(x,y). 

Then, because of the difference between the relations (14. ip and (14.21) we have in the right 
hand side of (14.61) an extra-term 2 5{x — y) A(x,y). One must collect all quasi-local terms 
appearing in this way and check if they can be put under the form cIqR(x, y) + i-£^R^(x, y) + 
■^pR^ty, x). If this can be done, then we can restore gauge invariance (at least for the tree 
contributions) by redefining the chronological products in an obvious way. 
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So the first problem is to find out the anomaly i.e. the expression appearing in the right 
hand side of (14.61) and the second problem is to see in which conditions it can be eliminated by 
a redefinition of the chronological products. Even the first problem is not exactly elementary 
in complex models as for instance the case of gravity; in [19] one can see for instance that not 
only terms of the type d^D^ix — y) A(x,y) can produce anomalies. So we need a systematic 
way to compute the anomaly. 

This suggests to make the following change in the description of the fields from section [U 
namely we replace the Pauli- Jordan distribution D m by some off-shell distribution which 
does not verify Klein-Gordon equation but converges in some limit (in the sense of distribution 
theory) to D m . For instance we can take 



where p m (A) is some function converging, say for A — > to the distribution 5(\ — m). In this 
way all the fields from Section [2] we become generalized free fields [16] i.e. they will verify all 
properties described there except Klein-Gordon and Dirac equations. For instance the off-shell 
scalar field we be denoted by $ , etc. However, for simplicity we will skip the index off if it is 
obvious from the context if we are considering on-shell or off-shell fields. 

If we keep the definion of the gauge charge unchanged we will loose the property Q 2 = 0. If 
we keep unchanged the expressions of the interaction Lagrangians from the preceding Section, 
but replace all fields by their off-shell counterparts, we also loose the relations (jl.fip . However, 
these relations will be replaced by 



with S 1 some polynomials which will be null in the on-shell limit. We will need these expressions 
in the following. In the following we will denote K c = K mc and we assume that all fields are 
off-shell (we do not append the index off). We have by direct computations the following results. 
A. The Yang-Mills Sector In this sector we have 

Theorem 4.1 The expressions S^ m have the following explicit form: 




(4.7) 



dnT 1 = i d u T^ + S 1 



(4.8) 



-0 = 
ym — 



i fabc U a V\ 




+ -u a u b K c u c j - i f' abc u a $ 6 K c § c 

t e a <g> Y 7 e - U a * M t e a <g> 7^ 7, * 

<8> Y 7, + u a ^I t~ e M ® Y 7e * 




-i u a <9„# 



(4.9) 



Also 




fabc U a U b K c v£ 



(4.10) 



and 



S 1 



0, |/| > 1. 



(4.11) 
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B. The Gravitational Sector In this case we have 
Theorem 4.2 The expressions S 1 have the following explicit form: 

S = = % re (2 u p dphap K h a ? + 2 dX Kp K h a ? - u» dph K h 

-<9X h K h + 2 d a u p h aP K h - 4 d p Ufl h up K W U + 2u p dpU u K u u ) (4.12) 

Also 

S» = -2i re u p d p u v K h pu (4.13) 

and 

S pu = i k (K u u u p d p u p + K u x dV u v ) - (/i <+ 1/) (4.14) 

C. Interaction between Yang-Mills and Gravitation In this case we have 
Theorem 4.3 The expressions S 1 have the following explicit form: 

S = S % = Ai f ab (w" F afJiU K b vl - Ai u, K a § a <%) + i f' cd u p d^ c K d $ d 

% — % — 

+- Up K^c e <g> 7 M 7 £ \1/ — — Up \I/c e ® 7 M 7e 



z - 1 

+- Up <W#c e ® 7^7 e * + - u M «9^M c £ ® 7 e * 

/ „..= « „„, 1 



+ - m m d'W ® 7„ 7, - - u M d^c~ e M ® 7, * 

z - 1 - 

-- u M cFW <g) 7^7^ <9^ - - n M c £ ® 7 e 

-~ Up ^c e <g> j v j e d u d p m + ^Up \I>c~ e M ® 7e d^tf (4.15) 

and 

= |/| > 1. (4.16) 

Now we proceed to higher orders of perturbation theory and we have the following result. 

Theorem 4.4 In the off- shell formalism we can choose the the second order chronological prod- 
ucts such that the following identity is true: 

d Q T(T I ^x),T I ^(y)) = iA-T(T^(x),T l2 (y)) + (-l)^^^,^) 

+T(S h (x),T^(y)) + (-1)^1 T(T^(x),S^(y)). (4.17) 

Indeed, if we make the substitution — > D^ ,off we obtain immediately the identity from 
(14.41) and (|4.8|) . Similar identities are true in the higher orders of perturbation theory. Let us 
consider the simplest case Ji = Ii = when we have 

d Q T(T(x),T(y)) = i-^T(T»(x),T(y)) + -^T(T(x),T»(y)) 

+T(S(x),T{y)) + {x^y). (4.18) 
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Now we have a very clear origin of the anomalies. It is elementary to prove that 

T(S(x), T(y)r cc = [ K m D^ s (x - y) A m (x, y) + d a K m D^ s (x - y) A^(x, y) 

m 

+d a d p K m D F m oS {x - y) A^{x,y) ] + ■■■ (4.19) 

where by • ■ ■ we mean terms where the Klein-Gordon operator (or Dirac operator) is acting on 
some off-shell field factor. So when we make the on-shell limit A — > we have 

T(S(x), T(y)f cc -> 6(x - y) A(x, y) + d a 6(x - y) A a (x, y) + d a d p 8{x - y) A^\x, y) (4.20) 

where the expressions A(x, y), A a (x, y), A^ a ^(x, y) are sums of the corresponding expressions 
A m (x, y), A^(x, y), Am^(x, y). In this way we have a systematic procedure to compute the tree 
anomalies in the second order of perturbation theory. For instance, the anomaly of the relation 
flUBD is 

A(x,y) = {5(x-y) A(x,y) + [dJ(x-y)} A a (x,y) + [d a dp5(x-y)] A^(x,y)}+(x y). (4.21) 

We investigate now in what conditions we can eliminate the anomaly by finite renormaliza- 
tions. The first trick is to use "partial integration" on the last terms with derivatives on the S 
distribution. We obtain the equivalent form: 



A(x, y) = 2 5(x - y) a(x, y) + 



_d_ 

dx c 



a a (x,y) + (x <r+ y) 



where 



and 



a{x,y) = A(x,y) 



_d_ 

dx~ c 



A a (x,y) + 



d 2 



dx a dx /S 



i a (x, y) = [d p 6{x - y)\ A<<*> (x, y) + 5{x - y) 







A a (x,y)--^A^(x,y) 



If we make the redefinition 

T(T»(x),T(y)) -> T(T»(x),T(y)) + i a"(z, y) 
of the chronological products we will put the anomaly in the form 

A(x, y) = 2 8(x — y) a(x, x) 

Now we have the following 

Lemma 4.5 The preceding anomaly can be eliminated iff the expression a(x) 
relative cocycle i.e. we have 

a = dnB - idmuB^ 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



a\x,x) is a 



(4.27) 
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for some Wick polynomials B and B^ L . The Wick polynomials B(x) and B^(x) are constrained 
by: (a) Lorentz invariance; (b) ghost number restrictions: 



gh(B) = 0, gh(B») = 1 (4.28) 
and (c) power counting which in our case gives: 

u(B) ,cu(B^) < 6. (4.29) 
The proof is very simple. Suppose that the anomaly (14.261) can be put in the form 

8{x-y) a(x) = d Q R(x,y)+i^R*{x,y) + -£jR' i (y,x). (4.30) 

with the expressions R(x, y), R^(x, y) quasi-local i.e. of the form 

R(x,y) = 8{x-y) B(x) + R^x,y) = 5{x - y) B"(x) + • • • (4.31) 

where ■ • • are terms with higher order derivatives on the 5 distribution. Then we immediately 
obtain from (I4.30p the identity from the lemma. Conversely, if the identity from the lemma is 
true then we take 



R(x,y) = S(x-y) B(x), i? M (x,y) = S(x - y) B^x) (4.32) 
and we have (I4.30p . 

So all we have to do it to compute the expression a(x,y) given by the formula (I4.23p . 
collapse the two variables to obtain the expression a(x) and impose the condition ( 14.271) . For 
simple models, as pure Yang-Mills theories, this computation is not very difficult but for more 
complicated models involving gravitation, the computation are very long and one can see the 
benefits of the off-shell method if one makes the comparison with the usual methods. 

In the same way one can treat the other identities of the type (14.171) i.e. for non-trivial sets 
of indices Ji, J 2 . 
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5 Second Order Gauge Invariance 

Now we turn to the question of gauge invariance of the model in the second order of perturbation 
theory. The case of Yang-M ills fields has been investigated previously [8] . One can eliminate 
the "tree" anomalies from the second order of the perturbation theory iff the following identities 
are true: 

f abcfdec fbdcfaec 4" fdacfbec 0, (5.1) 
fdcafceb ~ f dcbf cea = ~fabcfdeci (5-2) 
fcabfcde + fcdbfcae + fcebfcda = 0' V 5 G (5.3) 

[t e a ,tl\ = if abc t c , (5.4) 

t a S b ~~ S b^a = ^fbca S t ■ (5-5) 

For the case of pure gravity one can prove that the second order anomalies can be eliminated 
without further constraints [19]. Using the off-shell method one can compute the anomaly in 
a more systematic way. The anomaly can be eliminated without further constraints even in 
the case of massive gravity [12]. It remains to analyze the part of the interaction Lagrangian 
which describes the interaction between gravity and the other fields. The relevant expressions 
are T* nt , S( nt given previously. The main result is 

Theorem 5.1 The second order chronological products verify the gauge invariance condition 

d Q T(t int (x)Ant(y)) = i^r(CMMy)) + ^ T (W*),C(2/)) (5.6) 

in the second order of perturbation theory iff: 

fab = -- $ab, a,b G ii U/ 2 

fed = 2 5 cd} c,de h 

c £ = -i k (5.7) 

The finite renormalization of the chronological product T(t- mt (x),ti nt {y)) is given by the expres- 
sion R(x,y) = 2 S(x — y) B(x) where 

b =~\ k2 (V hT Tf F v „ + h 2 Fr F aprr - 8 V V K' K° + 8 K„ h F? F"„) 

+2 i J (y *„ *. - i H' 

-in f abc {-h v ail v bu F» u + 4v ap v bu \t p F v hp + 2u^ v» v v b d v u c 
-8m f bca (2 v% $ 6 d u <S> c + 2m b < v v b $ c + m b Ufl u a v£ $ c 

+™ fabc h $a $6 $ 

1 



-2i n 2 



4 



V h 2 ) ijj M <g> 7^ 
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+i u p u v (d p ip Ylei 1 + d v ip 7^7^ - i> l^led"^ - i> Yled^ip) 
+ l - u [H u p $ 7^7^ + % - {d p h va h a . p - d,u u u p ) i> 7^7^ 

~\ (V K + hh pu - Au p u v ) (d"$ Y 7^ - $ l v led^) 

+2i k < h pv $ t\ <g> Y 7^ (5.8) 

Proof: (i) We will compute the anomaly using the off-shell method described in the preceding 
Section. We first consider the contribution to the anomaly without Dirac fields. The expres- 
sion is very long (it has 75 terms!) and we observe that we can consider independently the 
contribution of the type x\ 22 Ui IJ2 and x y\ y-i 2/3 where Xj are gravitational variables and yj 
are Yang-Mills variables. So the expression a of the anomaly has two contributions of the two 
types described above. After a long computation one can put the first contribution in the form: 

a x = d Q B 1 + total divergence + 2 g ab W ab - 2 g' ab W' ab H (5.9) 

where the matrices g, g' are given by 

9 = 2/ 2 + « / 

g ' = 16 f + 8K / + (/ /)2_ 2 K/' (5.10) 
and the Wick polynomials W ab , W' ab by 

W ab = Au p d p u v F v a » d p u b - 8 dV h vp F pX F bX + 4 d p u p h p<7 F pX F b ° x 

+4 d p h pa F px F b ° x - u p d p h F bpry (5.11) 

and 

W' ab = u p (d p h^ d p <$> a d v $ b + hT d p $ a d p d u $ b ). (5.12) 

We have denoted by • • • in (15.91) the terms of canonical dimension less than 6 (these terms are 
monomials with mass factors) . The elimination of the anomaly (15.91) is possible iff we have 

9ab W ab - g' ab W' ab = d Q B[ + total divergence (5.13) 

From here we obtain 

9ab d Q W ab - g' ab d Q W' ab = total divergence (5.14) 

and if we make a generic ansatz for the divergence in the right hand side we can prove that we 
have g = 0, g' = i.e. 

2/ 2 + K / = 0, (f) 2 - 2k f = 0. (5.15) 

One can use the preceding relation to simplify considerably the forms of the matrices /, /'. 
Indeed, we already know that these matrices are real and symmetric, so they can by diagonalized 
with orthogonal matrices: 

f — C F C-\ f = a F' {C'r 1 (5.16) 
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where 

F ab = F a S ab , F' ab = F' a S ab . (5.17) 
From (I5.15P we have immediately 

2F 2 + k F = 0, (F') 2 — 2kF' = 0. (5.18) 

and if we substitute the preceding diagonal expressions we obtain 

F a = ~, F' a = 2n, Va (5.19) 

and we obtain the diagonal form from the statement. Moreover, if we compute all the terms of 
a\ (i.e. those with mass factors also) we obtain that we have 

ai = cLqB 1 + total divergence (5.20) 

with B\ given by the first two lines from the expression B from the statement, 
(ii) We turn now to the contribution a 2 and after a long computation we obtain 

a 2 = d Q B 2 + total divergence + ff « W abc - ~ + «£>) + & W$ + ■■■ (5.21) 



where the constants g are 



9al = fabe fed +(b^c) 
be = fed fdba ~ ^ fed fdba 

9 { l = 9 { l-^h d fL (5.22) 



and the Wick monomials W by 



= —An, V P J? 4. 7 ,P J? —Ah v F^ p F u 

VV abc — ^ U P U a r b r cpu T u p v a r b r cp°~ ^' t /if u a r b r c 



cp 



+h u a Fj; u F ctiV + 4^ u a F£ u d u u c (5.23) 

Wfl = 2 vP < d M $ 6 0„$ c - m m < 9„$ 6 0"$ c + 2 hT u a d^ b d u $ c (5.24) 

Wffi = Up, Fr $ b d^ c (5.25) 

and the dots ■ ■ • are terms of canonical dimension less than 6. However it is a remarkable 
fact that the expressions g^ bc , gf bc + g^ c and gf bc all null! This can be proved if we combine 
the diagonal form of the matrices /, /' with the other properties of the constants f abc , f' abc 
obtained in the analysis of the Yang- Mills model. Moreover if we compute the mass-dependent 
contribution in the expression a 2 we have after a long computation 

a 2 = d Q B 2 + total divergence + g abc $ a $ 6 <9 M $ C (5.26) 

where 

9abe = \ {-^fed + f' cd ) f'L- (5.27) 
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However, if we use as above the various properties of the constants /, /' and f" bc already 
obtained we can prove that the expression g a b c is completely symmetric so we easily show that 
the last term from (I5.26P is a relative coboundary. We obtain that a 2 has the form 

a 2 = dqB 2 + total divergence (5.28) 

with B 2 given by the lines 3,4 and 5 of the expression B from the statement. 

(iii) We consider now the Dirac terms of the anomaly. Again we have a very long list (73 
terms!) which can be group in two contributions ~ x± x 2 ip t/> and ~ x y ip ip where 
as above Xj are gravitational variables and yj are Yang-Mills variables. We start with the 
contribution 03 and after a long computation we can put it in the form 

a 3 = d Q B 3 + total divergence + 2 h pa A p{pr7} + 2«m /1 /i# (5.29) 

where 

+d p d a ip A e g> 7m lt i\) - i) A e ® 7m 7,^^-0 (5.30) 

A M = <9 M ?/> M A e ® j £ ip + ip M A e <g> 7e <9^ - $ A e <g> 7m le d p d a ^) (5.31) 
and we have defined the matrices 

A e = k d - i (c e ) 2 , e = ±. (5.32) 

We must have 

m m V A p{pa} + 2 /i A p = d Q B - id^ (5.33) 

and if we make an ansatz for the relative coboundary from the right hand side we immediately 
obtain that we must have A p ^ pa ^ =0, A p = and this leads to A e = 0. i.e. 

(c e ) 2 — — i k c e , Ve = ±. (5.34) 

We can use this relation to simplify the expressions for the matrices c e . Indeed, we have 
previously derived the self-adjointness property {c + Y = — c~ so the matrices 

c = l - (c+ + c-) d=\ (c+ - c-) (5.35) 

are verifying 

c f = -c, (c') f = d. (5.36) 

Because the matrix c is anti-self-adjoint we can diagonalize it with a unitary transformation: 
c = U C U^ 1 with Cjk = Sjk fj- Similarly we define d = U C U^ 1 so we have c e = U C e U" 1 . 
Now we have from relation (15.341) that 

(C 6 ) 2 = -1 k C e , Ve = ± ^ 
C 2 + (C) 2 = -i kC, C C' + C C = iKC' (5.37) 
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and this leads to the conclusion that the matrix (C) 2 is diagonal. Because all these matrices 
are finite dimensional, this means that the matrix C is also diagonal. If we substitute in ( I5.37P 
we get that the matrices C e are in fact proportional to the unit matrix. This fact is true then 
for the original matrices c e also. If we insert c e = X e Id, X e G C in (15.341) we get X e = —i k, Ve 
and this leads to the expression for c e from the statement. This condition ensures that 



where B% are the lines 6, 7 and 8 of the expression B from the statement of the theorem, 
(iv) Now we still have to analyze the contribution 04. After a long computation we get 



where B 4 is the last line of the expression B from the statement of the theorem. 

The expression for the counterterm B^ follows also from the preceding computations but 
we skip it. ■ 

6 Conclusions 

Second order gauge invariance does not lead to some new information concerning the Yang- 
Mills sector. We only get the diagonal structure of the matrices /, /' and c e for the interaction 
sector. This result can be interpreted as the universality of the interaction of gravity with other 
fields: one does not need new constraints to make the existence of the interaction Lagrangian 
possible. The extension of the method to the case of massive gravity is elementary and leads 
to the results from [T5] . 

In further publications we will consider loop contributions in the second order of the per- 
turbation theory and also higher order of perturbation theory using the same ideas. 



^3 = rfg-E>3 + total divergence 



(5.38) 



c-4 = dqB^ + total divergence 



(5.39) 
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